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Current-induced torques on ferromagnetic nanoparticles and on domain walls in ferromagnetic 
nanowires are normally understood in terms of transfer of conserved spin angular momentum be- 
tween spin-polarized currents and the magnetic condensate. In a series of recent articles we have 
discussed a microscopic picture of current-induced torques in which they are viewed as following from 
exchange fields produced by the misaligned spins of current carrying quasiparticles. This picture 
has the advantage that it can be applied to systems in which spin is not approximately conserved. 
More importantly, this point of view makes it clear that current-induced torques can also act on 
the order parameter of an antiferromagnetic metal, even though this quantity is not related to total 
spin. In this informal and intentionally provocative review we explain this picture and discuss its 
application to antiferromagnets. 



I. SPIN TRANSFER TORQUES 

The study of spin-transfer torques began in 1996 when John Slonczewskii and Luc Bcrgei— independently pre- 
dicted that magnetization dynamics can be induced by current in circuits containing noncoUinear magnetic elements. 
Berger's paper focused on the emission of spin waves as the source of magnetic excitations, while Slonczewski invoked 
conservation of spin angular momentum to infer magnetic dynamics. Slonczewski's observation that a net spin-current 
flux into a volume of magnetic material implies that a torque acts on the magnetization in that volume is the key idea 
for most theories of current-induced torques (CIT) . Over the past decade many experiments have confirmed^i'^'^'i'^'^ 
Slonczewski's predictions. There has also been theoretical progress, elaborating on Slonczewski's ideas and developing 
techniques which shed light on their implications for particular materials combinations and geometries. For example. 
Stiles and Zangwill explicitly exhibited all of the sources of net spin-current flux, namely spin-dependent transmission, 
spin precession, and rotation of reflected and transmitted spins Brataas et al. have formulated a general theory 
in which spin-dependent interface conductances, calculated from first principles or extracted from experimental data, 
can be combined using a generalized set of Kirchhoff laws to predict transport properties and magnetization dynamics 
in a circuit containing noncoUinear magnetic clemcntsiii Other approaches include solving generalized Boltzmam*i^ 
or spin diffusion equationsJ^ All of this theoretical work is directed towards evaluation of the net spin-current into a 
volume of material. Since the underlying systems possess conservation of total spin angular momentum, any theory for 
these systems must respect this global conservation law. For systems in which total spin is conserved, the relationship 
between net spin current and torques is very general. Appealing to spin conservation enables reliable predictions to 
be made about current-induced magnetization dynamics without having to specify which spins form the macroscopic 
magnetization or the microscopic origin of the effective magnetic fields which cause them to precess. The conservation 
laws assure that if one does the "bookkeeping" of spin properly, the current-induced torque acting on a volume may be 
inferred. We therefore sometimes refer to the point of view which utilizes powerful conservation-law consequences as 
the bookkeeping theory of spin-transfer. This point of view has had qualitative and quantitative success in describing 
experiments on spin-transfer in spin valves systems which are composed of several ferromagnetic nanoparticles, and 
on spin-transfer induced domain wall motion in ferromagnetic nanowires. 

The bookkeeping theory of spin-transfer torques raises two questions which often need not be answered explicitly 
and which initially drew our interest to the issue of current- induced torques. In attempting to provide answers to these 
questions that we find satisfying, we have been led to the theoretical picture of current-induced torques described 
below. 

i. What is the distinction in the spin-transfer picture between the electrons that carry current and the electron spins 
that compose the magnetic condensate? This question is particularly relevant for the transition metal systems which 
are the workhorses of metal spintronics since we know that both s-like and d-likc orbitals must be treated as itinerant. 
There is no clean distinction between the electrons which contribute the moments that order and the electrons that 
carry current. 

ii. Does spin-transfer occur in systems with strong spin-orbit coupling? Because of spin-orbit coupling, spin angular 
momentum is never really conserved, even in perfect crystals that have no disorder. How strong does spin-orbit cou- 
pling need to be to weaken or eliminate the spin-transfer effect? This question is particularly relevant to ferromagnetic 
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semiconductors like (Ga,Mn)As in which the spin-orbit coupling strength is^ comparable to the magnetic exchange 
energy and to the Fermi energy. 

II. CURRENT-INDUCED TORQUES AND NON-EQUILIBRIUM SPIN-DENSITIES 

In our picture spin-transfer torques arise as follows: The spin-density of electrons near the Fermi energy of a magnetic 
metal is altered when they carry a current through a non-coUinear magnetic environment. The change occurs as they 
realign their spins to sample the exchange field of the ferromagnet and thereby steer their spin orientations to match 
their non-coUinear environment. This change in spin-density leads to a change in the exchange field of the ferromagnet, 
which causes spins in orbitals far from the Fermi level to precess. We will refer to this as the current-induced torque 
picture of spin-transfer. Although consistent with the bookkeeping theory of spin-transfer for systems in which total 
spin is conserved, it suggests that the phenomena is more general. One implication as we discuss below is that 
current-induced torques drive order parameter dynamics in antiferromagnetic metals. 

A. Current-induced torques 

We first introduce some notation and provide a general orientation. The formalism we describe is applicable to any 
mean field theory, but in this article we consider primarily systems described by the local spin-density approximation 
(LSDA) of density functional theory (DFT). The notation we use below is appropriate for a Hamiltonian in a tight- 
binding or local orbital basis. We separate both the single-particle Hamiltonian and the density matrix into spin- 
dependent and spin-independent contributions: 

H.'. - Bf^^lA,,-r. (1) 

where f is the vector of Pauli spin matrices and i', i are site or orbital indices. In DFT, H^'i is the Kohn-Sham single- 
particle Hamiltonian. The notation for the spin-dependent part of the H is chosen to emphasize that it produces a 
spin-splitting A when it is orbital independent, as is often assumed in simple toy models of a ferromagnetic metal. 
(Note that m and A are in general complex for orbital off-diagonal elements, i.e. for i ^ i'; any Hermitian matrix 
may be uniquely written in the above form.) In LSDA, the interaction contribution to A and m are related locally at 
each point in space according to 

A(f) ~ AQ(n{f),rn{r)) rh. (2) 

where n and rh = mm are the local charge and spin-densities, respectively, and Aq is some parameterization of the 
exchange-correlation potential. Note that A(r) acts like an effective-magnetic field experienced by the Kohn-Sham 
quasiparticles. (In Eq. A(f) and m(r) are functions in real space, while in Eq. ([T|), Ai^i> represents the i, i' matrix 
element of the real space potential A(r), and ruiji the spin-dependent part of density matrix in orbital space.) A 
local relationship in space does not imply proportionality of the orbital representation matrix elements. In particular 
there are strong exchange interactions between s-like and d-like orbitals. 

We define the difference between spin-densities in the presence and in the absence of a current as the non-equilibrium 
spin-density, denoted by m*'' , and the corresponding difference in exchange-correlation potentials is denoted by A**" = 
Ao(?T.,m) m*'"/|m|. This expression assumes that \m*''\/m is small, something that is valid for any reasonable current 
strength. In a circuit with a noncoUinear magnetic configuration, the contribution to the local exchange-correlation 
effective magnetic field from the non-equilibrium quasiparticles will in general point in a different direction than the 
magnetic condensate, and this misalignment is responsible for the ensuing torque on the magnetic condensate - the 
spin-transfer torquei^^ The non-equilibrium spin-density driven by a source to drain bias voltage in a specific nano- 
scale circuit can be evaluated theoretically using non-equilibrium Greens function techniques, as we describe below. 
For a bulk magnetic metal with a smooth spin-texture, the non-equilibrium spin-density can be evaluated by treating 
both the spatial variation of magnetization direction and the electric field which drives a bulk current perturbatively. 

As an example, consider the spin valve structure of two ferromagnets separated by a spacer. If a bias is applied, 
there is a component of the non-equilibrium spin-density which is perpendicular to the plane spanned by the two 
ferromagnets. Because this spin-density is not parallel to the exchange-correlation potential A, the transport orbital 
spins precess as they move through the circuit to accommodate the change in exchange-field orientation. The mag- 
netization of both layers will precess around the local out-of-plane exchange field generated by this non-equilibrium 
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FIG. 1: Left panel: Ground state of a metallic ferromagnet. The low-energy collective degree of freedom, the magnetization 
direction, is the spin orientation of singly occupied orbitals. Right panel: Quasiparticles experience a strong exchange field 
A that brings majority and minority spins into equilibrium. Because this field is parallel to the magnetization it does not 
produce a torque. In an inhomogeneous ferromagnet, the spin orientation of the transport orbitals (in a window of width eV at 
the Fermi energy) must differ from the direction of the exchange-correlation potential A. The current-transfer torque is then 
produced by the transport-orbital contribution to the exchange field. 



spin-density, and this precession is the one which leads to current-induced magnetization switching (CIMS) (see Fig. 
©)• 

In the case where total spin is conserved, this view is essentially equivalent to the bookkeeping theory of spin- 
transfer, as can be seen by evaluating dm}''' /dt for a current-carrying quasiparticle with source to drain scattering 
boundary conditions: 

^ = ^ [m*^ H] = (Oi - 0„) + ixr [rc [m''' x a] ] . (3) 

The above relation can be obtained by evaluating the commutator directly and using the commutation properties of 
Pauli spin matrices. From Eq. ([T|), m*'' is the spin-dependent part of the scattering state density matrix and A is the 
spin-dependent part of the Hamiltonian, each of which is expanded in the 3 Cartesian components of Pauli matrices 
Tx,Ty,Tz. The trace in the second term is over orbitals in the subsystem (in the example shown by Fig. ([2|), the 
subsystem consists of all orbitals between planes 1 and n), and the spin-dependent part of matrices are multiplied in 
cross product form. 

The spin-current operator Qi is defined as the spin-current that flows between sites i and i + 1 , and is given by 

^ Tr [f {pkjHjk - HkjPjk)] ■ (4) 

k<i 
3>i 

In Eq. (j4|) the trace is over spin space, and Hij {pij) is the 2x2 spin matrix of the i, j-orbital component of H 
(p). The second term of Eq. ([3]) represents the torque present on the quasiparticle due to its misalignment with 
the magnetic condensate. In steady state transport, the left hand side of Eq. ^ vanishes, and so the divergence of 
the spin-current - or the net spin-current flux, is equal to the quasiparticle-condensate current-induced torque. The 
equation of motion for the magnetic condensate 'is^: 

— = -TO X A = -m*'' x A . (5) 
dt ^ ' 

From Eq. ([3]) , this implies that the torque on m can be calculated with either to*'' (the non-equilibrium spin-density) 
or (Qi — Qn) (the net spin-current flux), verifying the consistency between the approaches. In identifying the current- 
induced torques as arising from non-equilibrium spin-densities, we have accomplished more than simply rephrasing 
the bookkeeping argument. We have identiflcd the underlying microscopic interaction that is responsible for current- 
induced torques. The same mechanism is operational in systems in which spin is not conserved. In addition, this 
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which spin-currents Q are evaluated are between layers labelled and 1, and n — 1 and 



microscopic view allows for the evaluation of local torques on individual atoms which can also drive the order parameter 
of antifcrromagnets. The current-induced torque which acts on the order parameter of a volume of antiferromagnctic 
material is not related to the net spin-current flux into that volume. We comment more on these systems in Sec. IIVI 



B. Current induced torques and exchange interactions. 

So far we have identified the current-induced torque as resulting from the misalignment of the magnetic condensate 
with the exchange field contribution of non-equilibrium quasiparticles near the Fermi energy. This misalignment may 
seem out of place in view of Eq. ([2]) - in the ground state, the total spin and exchange field are aligned. However, 
even in equilibrium systems (systems that carry no charge current) with magnetic excitations (such as spin waves), 
the magnetic dynamics can be determined from an expression similar to Eq. ([5]) (assuming the excitation energy is 

"small" compared to other characteristic energies). In this case the Hamiltonian H' = - 1/2 (^A' ■ describing 

the magnetic excitation is constructed "by hand" starting from the ground state Hamiltonian (by, for example, 
applying a space-dependent spin rotation operator which describes a spin wave imposed on a coUinear ground state 

Hamiltonian). The resulting non self- consistent density matrix p' = 1/2 fp'"-' + m' • r ) can then easily be calculated. 



The torque on the magnetic system is then given by Tr Re to' x A' , as in Eq. ([5]), and yields proper values 

for magnetic exchange energies and spin wave dispersion relations. Indeed, it can be shown that this approach to 
calculating properties of magnetic excitations is equivalent to previous approachesii which calculate SE/6rh - the 
change in energy associated with magnetic excitations Srhr^ (As a technical note, we remark that the evaluation of 
6E/dTh, as well as our evaluation of torques, relies on an adiabatic condition for magnetic dynamics. The adiabatic 
approximation follows from the fact that electronic times scales are much faster than the time scales which characterize 
collective magnetic dynamics.—) Having established the close relationship between torques and variations in energy 
for excitations of equilibrium systems, one is naturally led to ask if such a relation exists for non-equilibrium systems. 
We address that point in Sec. IIIII 



C. How to calculate non-equilibrium spin-densities 



We now briefly describe the technique we use to evaluate spin-densities and provide some references. This in some 
sense represents a technical point outside the focus of this article (although certainly a point of practical importance) . 
We adopt a Landauer-Buttiker approach for describing non-equilibrium transport. In this approach, the bias voltage 
is represented by placing the system in contact with particle reservoirs with chemical potentials fj.s = cf + eVs/^ 
and = — eVB/2 in source and drain respectively. When the system Hamiltonian and its coupling to source and 
drain electrodes is time- independent (an assumption which follows from the adiabatic approximation described above), 
electrons with energies inside of the transport window fijj < E < fig solve a time- independent Schrodinger equation 
with incident-from-source scattering boundary conditions. In the Landauer-Buttiker approach, the conductance is 
simply proportional to the quantum-mechanical transmission probability of the scattering state. There are a number of 
techniques available to solve this system. We employ a non-equilibrium Green's function (NEGF) methodic combined 
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with density functional theoryi^ The NEGF formalism yields all quantities of interest: the conductance, and the 
contribution to the density matrix from the non-equilibrium scattering states, from which m*'' is determined. A useful 
introduction to the NEGF formalism can be found in Ref. Ull, and a more formal account is given in Ref. [l^ • More 
details of how NEGF can be used to find current-induced torques can be found in Ref. [2^. 



III. REACTIVE AND DISSIPATIVE TORQUES 



In the previous section we described how current-induced torques arise as a consequence of the interaction between 
the magnetic condensate and the spin-density of non-equilibrium quasiparticles, a point of view that is an extension 
of the notion of spin-torques arising from spin-currents. The exchange-correlation field of a non self-consistent spin- 
density may also be associated with the torques arising from magnetic stiffness, and in this context these torques 
are related to the variation of an energy functional. We now comment on the possibility of finding current-induced 
torques from an energy functional, which is of particular relevance in view of the ongoing discussion regarding the 
form of magnetization dampingi^i^ 

The Landau-Lifshitz- Gilbert equation describes the dynamics of the magnetization. When spin-transfer torques 
are included, it is: 



dm 
'dt 



SE[m] 
hSrh 



am X 



dm 

'at 



dm 
'dt 



(6) 



STT 



Here rh is the magnetization vector normalized to unit magnitude. The first term on the right hand side describes 
precession of the magnetization around an effective field, written here as the functional derivative^ of the energy. In 
the so-called micromagnetic theory used to describe smooth spin- textures, the energy functional contains contributions 
from Zeeman coupling of the magnetization to an external magnetic field, anisotropy energy contributions due to 
spin-orbit coupling and magnetostatic interactions which violate spin conservation, as well as the energy cost of 
magnetization variation which is referred to in the context of micromagnetic theory as the exchange energy. (As 
experts are well aware, the theory of magnetism refers to many different things as exchange energies. The many uses 
of this word is fitting, since magnetism is always intimately related to electronic exchange processes, but it can be 
confusing.) In Eq.([6l), the term proportional to a is the Gilbert damping term. 

For smooth spin-tcxturcs, the spin- transfer torque terms can be expressed in terms of the leading order of a gradient 
(spatial derivative) expansion: 



dm 
'dt 



-(vs • V)?7i + (3m X (vg • V)to , 



(7) 



STT 



where the velocity Vs = —o?P]/\e\ is proportional to the electric current density j. Here a is introduced as the lattice 
constant of a fictional lattice of unit magnetic moments representing a magnetization density a^^ , and P is the spin 
polarization fraction of the current. Finally, the electron charge is denoted by — |e|. The first term in Eq. ([7]) is known 
as the adiabatic spin-torquc. When v,, is defined as above, this term is just the gradient expansion limit of the spin- 
conserving spin-transfer torque discussed in the previous section. The second term proportional to the dimcnsionlcss 
parameter (3 is commonly (and perhaps inappropriately - see below) referred to as the nonadiabatic spin-torque. Since 
m is a unit vector, it is perpendicular to its time derivative m. It follows that when both P and (3 are treated as 
phenomenological parameters, Eq. ([7]) is quite general, assuming only that the current-induced torque is linear in 
current. From a symmetry point of view, the first spatial derivative spin-transfer torque terms in the generalized 
Landau-Lifshitz equations are allowed because current-flow breaks inversion symmetry. The micromagnetic exchange 
term, which is proportional to the second spatial derivative of the magnetization, is the leading order term in a 
gradient expansion in the absence of current. Eq. ([7]) can be derived microscopically,i2Si^i22ii2^ including the spin- 
transfer torque terms, by using non- linear response theory to describe the response of magnetization to external fields 
in the presence of a transport current. 

The analogous equation for the free nanomagnet in a spin- valve system is 



dm 
'dt 



-gjm X (rhp x m) — f3gj [rhp x m) 



(8) 



STT 



with rhp the magnetization direction of the pinned magnet. In this equation g ~ a!^Pl{2\e\t) with t the length of the 
free ferromagnet in the direction of current flow. In this section we follow one common usage by referring to the first 
term as the Slonczewski spin-torque, and to the second term as the effective-field spin-torque; physically g represents 
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the efficiency of the Slonczewski spin-torquc. An important observation is that Eqs. ^ and ([8]) are two sides of the 
same coin. If we put rhp = Th{r — dr) in Eq. ([5|), with r the coordinate in the direction of the current, we obtain, after 
expanding to lowest order in dr the result in Eq. ([7|). In the case of spin- valve structures we know from microscopic 
theory that the approximate expression for g has corrections from spin-dependent reflection and other effects that are 
easily captured by microscopic interface calculationsJ^iii 

A useful classification of the terms in Eq. ([6|) follows from an examination of how they behave under time reversal 
operations. Letting t i-^ —t, we have that rh,rhp i-^ —rh, —rhp, and Vs,j i—^ — v^,— j. Furthermore, E[—m] = E[rn\. 
After implementing this operation in Eq. ^ we observe that there are two kinds of terms, reactive terms that are 
even under time reversal and dissipative terms that are odd. (Readers should be aware that there exist different points 
of view regarding this issue, found in Ref. HO and [l^.) The terms reactive and dissipative originate from linear 
response theory in which they play a similar role, distinguishing response that is in phase with a periodic driving 
force from response that is out of phase and therefore dissipative. The first and second terms on the right hand 
side of Eqs. ([7]) and ([8]) respectively represent corrections to the reactive and dissipative terms in Eq. ([6]) due to 
current flow. In the smooth texture (continuum) limit this is obvious from a microscopic point of view, because the 
terms proportional to a and /3 emerge from microscopic linear-response theory as dissipative parts of the spin-density 
spin-density response function in the presence of current j2i To many readers the classification of the Slonczewski 
spin-torque as reactive may come as a surprise: after all, the Slonczewski spin-torque competes with the Gilbert 
damping in current-driven magnetization reversal. The latter phenomenon is dynamic, however, and does not imply 
that the Slonczewski spin-torque is dissipative. This classification is consistent with our picture that current-induced 
torque phenomena primarily reflect precession around a transport-induced contribution to the exchange-correlation 
effective magnetic field. 

The fact that both the adiabatic spin-torquc and the Slonczewski spin-torque are both reactive triggers the question 
whether they can be derived from an energy functional. Indeed, as we show below, the action 



A[m] 



dt 



hVs^-^A^i {rh(x, t))V^m^i (x, t) + hA{m{x., t)) 



drh{x, t) 



dt 



E\m] 



(9) 



with A the vector potential of a magnetic monopolc determined by 



drhyl 



(10) 



that enforces spin quantization,— reproduces, upon variation, the equation of motion in Eq. ([7|) for a = (3 = Oj^ 
We note that in the above equation e^^'^" is the antisymmetric Levi-Civita tensor and sums over repeated Cartesian 
indices 7, 7', 7" G {x, y, z} are implied. 

To obtain the equation of motion corresponding to the action in Eq. ([9]) wc have to calculate the variation of the 
action, i.e., SA/Sm^ and put it equal to zero. We find in first instance that 







SA _^/dA^ 



Sm 



drh^ 



d_ 

dt 



d 



SE 
5rh^ 



(11) 



Using Eq. (fTO|) we have that 



dA, 



^ _ dAy 



(12) 



which we use to rewrite Eq. pT|) as 



hrh X 



dt 



+ (v, • V) 



SE_ 

dm 



(13) 



Taking the cross product of the above equation with m and using the fact that tti is a unit vector we obtain 
Eqs. dSD and dZl) for a = /? = 0. 

The term in the action in Eq. ^ that is proportional to Vg is physically understood as the Berry phase acquired 
by the electrons as they drift through a non-coUinear magnetization texture. Such Berry phases occur naturally in 
spin systems. To see this in more detail consider first the simple case of a spin S in a time-dependent unit magnetic 
field m{t) with hamiltonian Ti. = —rh{t) ■ S. Suppose the magnetic field is varied very slowly from m(ti) to m{tf), 
with rh(tf) = rhiti), so that the system remains in its ground state. One can show^l that the quantum mechanical 
wave function of the spin acquires a nontrivial phase factor e~*^^. Here, f2 is the area on the unit sphere enclosed by 
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the path that the spin traces out as it adiabatically follows the magnetic field. Using Stokes' theorem we write this 
area as a line integral over the boundary of fl, that is, 



/ m-da= / [V,h X A(m)] ■ da ^ A(m) ■ dl , (14) 
Jn Jn J on 

with the monopole vector potential A(?7i) determined by Eq. (fTUl) . Note that the above also shows that although the 
vector potential is only defined up to a gauge transformation A ^ A — VmA, with A an arbitrary function of m, the 
physical quantity of interest, namely the area fi, is unaffected and well-defined. Using the above results we observe 
that the term proportional to in the action in Eq. ([9]) is determined by the area that the magnetization traces out 
on the unit sphere. 

Having found an action that reproduces the equation of motion for the magnetization including the adiabatic spin 
transfer torque, we define an energy functional 

f dx ^ 

Ej[m] = / —hvs^^A^r{m{x,t))W^m^'{x,t)+E[rh]. (15) 



Using this energy functional the equation of motion for the magnetization direction is written as 

dm ^ f SEj ' 



dt V hSifi 



(16) 



To illustrate that the energy functional in Eq. psp is indeed the energy that is minimized in the presence of current 
and as such a useful concept, we now add dissipative terms to the above equation. In particular, we consider the case 
a ^ and f3 = 0. In that case we have that 

dm ^ [ SEj \m] \ ^ dm , „, 

so that, since a > 0, the energy functional Ej[m] is indeed minimized according to the above equation of motion. 
This energy functional is relevant for understanding the so-called "intrinsic pinning" of a domain wall that occurs 
for P ~ 0. "Intrinsic pinning" refers to the fact that a magnetic domain wall in a perfect material is displaced only 
for a sufficiently large applied current, and is "pinned" otherwisei^ Minimization of this energy functional yields a 
physically clear explanation of intrinsic domain wall pinning i^i^l 

We end this section with some comments regarding spin valves. Writing down an action that reproduces Slon- 
czewski's spin-torque directly for the spin valve case turns out to be more complicated because of the "discretization 
limit" one has to take in going from a smooth magnetization texture to a spin valve. We speculate, nonetheless, that 
the above action in principle is sufficient, as it reproduces Eq. ([7]) and we have seen the equivalence of Eq. ([7]) and 
Eq. ([5]). At this point some readers may object that it is in fact the dissipative effective field torque proportional to 
f3 in Eq. ([5]) which can be derived from the energy functional 

E[m] ^ -Pgjrhp ■ m . (18) 

We note, however, that this energy is odd under time reversal and therefore not a proper energy. Only when the 
pinned magnet is kept fixed as t i— > — < it is even under time reversal . This corresponds to regarding the pinned 
magnet as distinct from the dynamical system at hand (the free magnet). 

Among the many questions which remain we mention two. i) Is the energy Junctional from which the reactive 
spin-torques are derived a practical concept? The prefactor of the monopole vector potential contains the spin-current 
evaluated in the coUinear situation. Is this useful for calculations of the Slonczewski spin-torque efficiency? We note 
that a different perspective on this issue can be found in Ref. Issl . ii) Can useful predictive expressions he derived for 
the dissipative coefficients a and P of real ferromagnetic materials ? 



IV. CURRENT-INDUCED TORQUES IN FERRIMAGNETS 



For systems in which spin is approximately conserved, such as transition metals, the current-induced torque picture is 
complementary to the standard bookkeeping argument of spin-transfer, as described earlier. For systems in which spin 
is not conserved (i.e. systems with strong spin-orbit coupling), the bookkeeping argument is no longer valid, and the 
current-induced torque picture must be adopted. Examples of such systems include diluted magnetic semiconductors 
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(DMS) and rare-earth elements. The current-induced torque picture has been previously employed for these systems, 
for example by Nguyen et al. who calculate the current-driven domain wall mobility in the DMS (Ga,Mn)Asi^ 

Recent experiments on ferrimagnets provide an interesting test of current-induced torque theories. Jiang et al. 
consider magnetoresistance and CIMS in a system consisting of a ferromagnetic CoFe fixed layer and a ferrimagnetic 
CoGd free layeri^ In CoGd, the magnetization of the two sublattices Co and Gd point in opposite directions. For 
temperatures above the magnetic compensation temperature Tmc, the net magnetization points in the same direction 
as the Co sublattice, while for temperatures below TmCi it points in the direction of the Gd sublattice. The net 
magnetization vanishes at T = Tmc- 

Transport is dominated by orbitals from the Co sublattice, so that the magnetoresistance probes the relative 
orientation of the Co sublattice with the ferromagnetic pinned layer. Therefore the magnetoresistance changes sign 
at T = Tmc (so that for T < TmCi the magnetoresistance is negative - in this regime the Co sublattice is antiparallel 
with the net magnetization). 

The experimental result of key interest is that the sense of the CIMS changes sign not at T = Tmc, but at a 
higher temperature. The interpretation provided for this experimental result was that at all temperatures, the CIMS 
is determined by the relative orientation of the total angular momentum of pinned layer and free layers. The angular 
momentum L and magnetization M of a sublattice are related hy L = M where ^ = g^s/fi is the gyromagnetic 
ratio, whose g-value depends on the the degree of spin-orbit coupling. Making the assumption that different (^-factors 
can be assigned to the Co and Gd sublattices, we can see that CoGd has the curious property that in a certain 
temperature range, its magnetization M = Mco + Mq^ and angular momentum L = A/co/tco + MQd/jod can have 
opposite signs. 

At temperatures T > Tmc, the magnetization and angular momentum are parallel, and are dictated by the Co 
sublattice, while at very low temperatures, the magnetization and angular momentum are also parallel and dictated 
by the Gd sublattice. In these regimes, let us suppose that a positive current leads to parallel alignment of the 
CoGd magnetization with the fixed layer magnetization. In the intermediate temperature range where M and L have 
opposite sign, so the interpretation goes, a positive current aligns the angular momentum L (not M) with the fixed 
layer magnetization. 

We believe that this is an important experimental result for understanding the basic nature of current-induced 
torques. The interpretation provided by the experimentalists extrapolates the bookkeeping theory from one based on 
conservation of only spin-angular momentum to conservation of total angular momentum. We do not believe that 
this extrapolation is well motivated. For one thing, orbital angular momentum is not even approximately conserved 
in a crystal. We do not believe that the temperature at which the sense of CIMS changes sign must equal the 
temperature at which the total angular momentum changes sign. It is true that because of spin-orbit interactions, 
ferromagnets do have an orbital angular momentum, but neither spin nor orbital angular momentum are conserved. 
Indeed, understanding the magnitude of the orbital contribution to the magnetization of metallic ferromagnets has 
proved to be a challenging problem for many-body theoryi^ 

Our microscopic picture of current-induced torques provides some guidance on how to think about these inter- 
esting experiments. A systematic inclusion of orbital magnetism is provided by current-density functional theory 
(CDFT)i^i^ CDFT is an extension of DFT in which the current density (as opposed to charge density) plays the 
central role, and is able to treat many-body systems in magnetic fields of arbitrary strengthj^ The exchange-correlation 
potential is more complicated in CDFT than in DFT, and includes a vector potential. However, a standard practice 
is to simply employ density functional theory with the standard exchange-correlation potentials that depend only on 
spin-density, even when spin-orbit coupling is included in the Kohn-Sham single-particle equations. One formal justi- 
fication for this comes from the relativistic spin-density-functional formalism4i This formalism allows for relativistic 
corrections to the exchange-correlation potential, but the corrections do not have an overwhelming effect. Taking 
this as a starting point, the current-induced spin-torques in CoGd would be proportional to induced spin-densities 
and would be influenced by the spin-orbit interaction terms in the Kohn-Sham equations for the current-carrying 
quasiparticles. 

Another aspect that appears to play a crucially important role is the thermal fluctuations of the magnetization, 
something which has not usually been accounted for by theory. The effect of these fluctuations is not clear a priori 
(and indeed represents an interesting avenue of research in its own right). At finite temperature, both the Gd moments 
and the Co moments will fluctuate in orientation, so that the magnetization cannot be assumed to be coUinear even 
within the CoGd nanomagnet. Presumably the change in sign of the average magnetization occurs for T > Ta/c 
because the Gd moment orientations fluctuate more strongly. 

It is interesting to consider how we would expect thermal fluctuations in such a system to influence macroscopic 
current-induced torques. Is the explanation for these experimental results related only to thermal fluctuations, with 
spin-orbit coupling playing an inessential role? Since the sum of the current-induced spin-torques on individual atoms 
will not in general be perpendicular to the total spin-density, its effectiveness in driving spatially coherent precession 
of the typical non-coUinear spin-density is not simply related to the net spin-current through the nanomagnet. The 
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intriguing experimental results in CoGd may indeed be pointing to a non-trivial role for thermal fluctuations of the 
spin-density in current-induced torque phenomena, when these fluctuations are large. 

V. SPINTRONICS IN ANTIFERROMAGNETS 

The more general nature of the current-induced torque picture suggests that it should be operational in more general 
circumstances, for examples in materials with more complex magnetic order than simple coUinear ferromagnetism. We 
have recently considered^'^i^'^i^ the effect of current-induced torques in circuits containing different combinations 
of antifcrromagnctic (AFM) and ferromagnetic materials. To date, the role of antiferromagnets in spintronic devices is 
to pin a ferromagnetic layer's orientation via an effect known as "exchange bias" We propose that antifcrromag- 
netic materials can serve as building block for circuits which display effects such as giant magnetoresistance (GMR) 
and CIMS, with the staggered order parameter of the AFM playing the role of the orientation of the FM. We briefly 
discuss our results for the sake of illustrating the qualitatively new features than can emerge from antiferromagnetic 
spintronics. These early efforts may help point the way to fruitful directions to consider in moving forward. 

The property that the local spin-transfer torque on each atomic site is equal to the net spin-current into that site 
is equally true for ferromagnets and antiferromagnets. The difference between the two-cases is the nature of the 
magnetic order. The magnitude of current-induced torques is nearly always dwarfed by the size of the torques caused 
by equilibrium exchange interactions when the relative orientations of different moments is distorted. In the case of 
a ferromagnet, equilibrium exchange orientations keep all the moments in a nanomagnet essentially rigidly parallel; 
the current-induced torques are generally not strong enough to change these relative orientations. The sum of all the 
current-induced torques acts on the overall orientation of all the moments. For this mode of magnetic dynamics the 
current-induced torque only has to compete with much weaker anisotropy, magnetostatic, external field, and damping 
effects. In an antiferromagnet, the mode of magnetic dynamics on which the current-induced torque can have an 
effect is the rigid motion of all the moments of the antiferromagnet. For an antiferromagnet with two sublattices with 
opposite orientations, it is the difference between the sums of the torques over the two sublattices which has an effect. 

GMR and current-induced torque (GIT) effects in ferromagnets rely on the interplay between electron transport 
and magnetic order. The source of this interplay is the strong spin-dependent exchange-correlation potentials seen by 
current-carrying quasiparticles, which is the result of a spin-dependent Fermi surface. Antiferromagnets do not posses 
a spin-dependent Fermi surface, the characteristic that is so essential to conventional spintronics, so it is clear that 
any GMR and GIT effects in AFM circuits must have a fundamentally different origin. (Indeed the only qualitative 
imprint of the ordered state on the electronic structure of AFM is the formation of a gap at the spin-density wave 
vector; we denote the spin-density wave vector by Q .) We can classify our studies by the relative orientation of the 
current J and the spin-density wave vector Q. 

Case 1: J \\ Q. In Ref. E we considered an antiferromagnets in which the exchange splitting changed sign on 
alternate lattice sites in the current-motion direction (J || Q). The structure we studied with this kind of model was 
similar to that of a normal spin valve, except that the magnetic nanoparticles separated by a normal metal spacer 
were antiferromagnetic rather than ferromagnetic. We refer to this type of structure as an antiferromagnetic spin 
valve. The model we studied in this paper was a single-band toy model. In the J || Q case the interface layer of 
the AFM is uncompensated, i.e. it has nonvanishing total spin. In this model, we find a difference in conductivity 
according to whether or not the AFM layers adjacent to the spacer are parallel or anti-parallel, an effect we refer to as 
antiferromagnetic giant magnetoresistance (AGMR). We also find the remarkable property that when the AFM layers 
are noncoUinear, the out-of-plane spin-density in the AFM is exactly constant in our lattice model antiferromagnet 
and exactly periodic in a continuum model antiferromagnet. Since the out-of-planc spin-density is responsible for 
the current-induced torque, the torque is also constant in magnitude and alternating in sign throughout the volume 
of the AFM. By itself, this property implies that the critical current for switching is independent of the AFM layer 
thickness. This is in stark contrast to the FM case, where the spin-torque decays rapidly away from the interface as 
the result of interference between different transverse channels' spin-density. 

Since the staggered magnetization of antiferromagnets is not conserved, the staggered torque which drives order 
parameter dynamics is not protected by robust conservation laws. Indeed, both the AGMR and the current-induced 
torques of the uncompensated AFM spin-valve toy model can be seen to follow from phase-coherence effects: specif- 
ically the difference in the phase acquired by up and down spins as they traverse the circuit in various transverse 
channels. The presence of the current-induced torque in this geometry may also be understood from the bookkeeping 
perspective: the difference in phase acquired by the up and down-component of an electron spin as it traverses a 
single uncompensated layer results in a net spin flux into that layer. 

The toy model calculations demonstrate that GMR can occur in purely antiferromagnetic spin valves and that 
current-induced torques can drive antiferromagnetic order parameters. At the same time, the toy model considerations 
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also suggest that in the AFM case the effects are more easily weakened by inelastic scattering and more sensitive to 
the details of the various interfaces in these layered systems4^ 

The toy model does not fully capture all the physics that is present in realistic AFM spin-valve structures. To 
partially assess the significance of real-world complications we have performed ab initio calculations for a system 
consisting of antiferromagnetic Cr leads separated by a Au spacer 4i The AGMR in this calculation is not primarily 
due to phase-coherent effects, but rather to spin-polarized interface resonance. The [001] surface of Cr with spin 
density wave in the [001] direction exhibits very different properties than that of the bulk, with an enhanced magnetic 
moment and a spin-dependent local density of states at the surfacc4^ The spin-dependence of the surface density of 
states at the Fermi energy results in a spin-polarized current as electron flow passes through the surface. The presence 
of Au spacer has little effect on this surface state, and the Cr layer can essentially be thought of similarly as a FM 
layer, with the surface magnetization playing the role of the FM magnetization. 

Another different aspect of this system is that the non-equilibrium out-of-plane spin-density is not periodic as in 
the simpler models, but partially decays away from the interface, because of the complex Fermi surface of Cr. For 
many transverse channels, the spin-dependent scattering state is a linear combination of Bloch states with different 
wave- vectors in the transport direction kz . The spin-densities of these scattering states then show an oscillatory spatial 
structure, and averaging over the Fermi surface results in destructive interference of the non-equilibrium spin-density 
away from the interface, as in the conventional FM case. We address some open questions regarding this system at 
the end of this section. 




(a) (b) (c) 

FIG. 3: Current-induced torques due to a compensated antiferromagnet. The arrows above the structure indicate the electron 
flux spin direction. The white arrows indicate the ensuing current-induced torques on the FM. The torque vanishes at ^ = 90° , 
and varies as sin 20. 

Case 2: J J- Q. We have also considered a system consisting of a ferromagnet and antiferromagnet in which 
the current is perpendicular to the spin-density wave ( J _L Q). The total magnetization of each layer in the 
current-carrying direction vanishes in this case, i.e. the magnetization is compensated. Most AFM materials used 
in magnetoelectronics are fully compensated, i.e. the spin-density sums to zero in every lattice plane perpendicular 
to the current direction, or at least nearly so. Direct interfaces between nearly compensated antiferromagnets, which 
perform the exchange bias function, and ferromagnets are common in spintronic circuits. The current-induced torques 
discussed in this section can be dominant when the antiferromagnetic and ferromagnetic layers are separated by a 
spacer which reduces the importance of direct exchange interactions between the two actors. Symmetry considerations 
imply that the current-induced torques between ferromagnets and compensated antiferromagnets differ qualitatively 
from the torques between ferromagnets. 

The total current-induced torque acting on a FM nanoparticle can alwaya^^ be expressed in terms of the difference 
between incoming and outgoing spin-currents. The presence of a ferromagnet will in general induce a nonzero spin- 
current at the AFM-FM interface. When spin-polarized electron flux from the AFM with orientation uafm enters 
a FM with orientation fipMj the spin-current entering the FM will have some component in the fijipM direction. 
It follows that, just as in the familiar case where both materials are FMs, a current-induced torque will act in the 
plane dcflned by hafm and fipM, as illustrated in Fig. ([3]). (Out of plane torques are also non-zero but tend to 
be much smaller.) Spin rotational invariance of the overall circuit implies that the in-plane torque must be an odd 
function of the angle 9 between hpM and uafm and that it can therefore be expanded in terms of a sine-only Fourier 
series, vanishing for both parallel and antiparallel coUinear conflgurations. In this case, reversal of the AFM moment 
direction is equivalent to a lateral translation which cannot influence the current-induced torque. It follows that 
in the compensated AFM case the torque is invariant under 9 — > 9 + tt, restricting its Fourier expansion to terms 
proportional to sin(2??6'). It follows that the torque vanishes when fipM is perpendicular to tiafm, and undergoes 
a sign change for 6* n — 9, as illustrated in Fig.( [S]). The property that the torque acting on a FM due to a 
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compensated AFM vanishes not only for coUinear but also for perpendicular orientations is primarily responsible for 
a novel current-induced-torque phase diagram, which differs drastically from the now familiar applied- field/current 
phase diagram for FM layers.— 

We emphasize that in making this general argument of the symmetry-constrained form of the current-induced 
torque on the FM, we have appealed to the bookkeeping picture. The bookkeeping picture combined with the 
assumption that the spin current becomes aligned to the local magnetization represents a simplification of great 
conceptual and practical utility. The calculation of the total current-induced torque in a spin valve then requires 
only the determination of the spin-current in the spacer. On the other hand, for the AFM it is necessary to find the 
staggered torque, and this convenient picture of finding net torques from a single spin-current value does not apply. 
To date, we have found torques on AFMs by adding up torques on individual atoms, and their form and magnitude 
seemingly can not be so easily be anticipated a priori. 

We have performed a realistic calculation of the torques present between when ferromagnetic Co is adjacent to the 
antiferromagnetic compound NiMn. As expected on the basis of the symmetry considerations explained above, the 
current-induced torque is of sin 26 form. We find that the current-induced torque efficiency acting on both layers is 
substantial and of the same magnitude as that found in common FM systems. For electron flow from FM to AFM, 
the CIT tends to align the axis of the AFM with that of the FM, and to make the FM perpendicular to the axis of 
the AFM. Conversely, for electron flow from AFM to FM, the CIT tends to align the FM with the AFM axis, and 
make the AFM axis perpendicular to the FM (within their common plane). Put another way: the current-induced 
torque tends to drive the orientation of the downstream material (AFM of FM) parallel with that of the upstream, 
and to drive the upstream material orientation perpendicular to the downstream. 

The effect of such a torque on the stability of the FM orientation is most unusual for electron flow from FM to 
AFM. If the FM is an easy plane ferromagnet (hard axis = x), and the AFM axis (AFM axis = z) is assumed to lie in 
the ferromagnet's easy plane and be fixed, then for sufficiently large current, the stable configuration for the FM is to 
point approximately out of the easy plane. Fig. ([4]) shows the stable magnetization orientation phase diagram versus 
applied field and current. The applied field is scaled by the demagnetization field of the FM, h = Happ/ Hdemag, and 
the conversion of the dimensionless current-induced torque hci into a real current J (assuming a demag field of 1 T) 
is J = [he it) X 3.8 • lO^A/cm^, where t is the thickness of the FM layer in nm; hci < represents particle flow from 
FM to AFM. 




-2-1012 -2-1012 
h h 



FIG. 4: Magnetic configuration (Mx,Mz) and peak of power spectrum Pz (arbitrary units) versus applied field and current. 
Also shown is stability boundaries found analytically (the labels ±x, ±z refer also to solutions which point approximately in 
these directions). The stability boundary plot also shows the reduced out-of-plane solution space for negative to positive field 
sweep with a dashed line. 



A stack design which may be employed in future studies to investigate this effect is shown in Fig. The AFM 
should be pinned in some way, here we've indicated pinning with the exchange bias effect by placing an larger FM 
adjacent to the AFM. As discussed in Ref. H^, it should not be necessary for the AFM to be single domain in order to 
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see the effect. The orientation of the free FM may be detected via the GMR effect with the pinning FM. The pinning 
FM may play a role in the stability of the free FM, but its role is sufficiently distinct from that of the AFM that 
it may be possible to nevertheless detect the AFM current-induced torque. A detailed discussion of this unexpected 
behavior can be found in Ref. IZsl . The qualitatively distinct nature of the applied current-applied field phase diagram 
suggests that there may be other interesting aspects of this form of torque to be explored as well. 



Applied field 
► 

Free FM 
N 

AFM 



Pinning 
FM 



FIG. 5: Possible stack configuration to study tfie influence of compensated AFM on FM orientation as a function of applied 
field and current. The AFM is assumed to be fixed by the larger pinning FM via the exchange bias effect. 

What next for AFM spintronics? There are other permutations of J, Q orientation and FM/AFM stack design to 
be explored. We believe the key step necessary for making further progress is finding materials systems which clearly 
demonstrates the effect of current-induced torques on or by an AFM layer in a way which can be distinguished from 
FM STT. 

Experimental challenges abound for studying AFM spintronics. Chief among them is that the structure and 
orientation of an AFM is difficult to measure and difficult to control. The easiest way to control the AFM orientation is 
with an adjacent FM via the exchange bias effect. However this FM layer may exert spin-torques of its own, potentially 
obscuring the role of the AFM. It would be preferable therefore to avoid using a pinning FM, or to somehow remove its 
effect on transport. To the extent that effects rely on flat interfaces or phase-coherence, the experimental challenges 
become more severe. Nevertheless, recent experiments^ established a dependence of unidirectional exchange bias 
fields on current, providing indirect evidence that current-induced torques are present in AFMs. Other experiments 
attribute observed changes in exchange bias, steps in differential resistance, and the statistics of thermally activated 
switching to current-induced torques on the antiferromagnct of an exchange-biased spin valvcj^ 

The role of disorder on both toy and more realistic models is likely much more important in antiferromagnets 
than in fcrromagncts. This role has not yet been realistically assessed, and represents an important step forward in 
determining the viability of AFM spintronics which can be achieved by theoretical work alone. For example, some 
important properties of the J \\ Q geometry seem to rely to some extent on a clean, uncompensated AFM surface. 
The robustness of the GMR and CIT when this theoretical assumption is relaxed is still largely unknown. 

Finally, so far we've only considered the effect of current on AFMs in multilayer geometries. It remains to study 
the effect of current on continuous AFM textures and domain walls, although some work is currently underwayi^ 
Domain walls in AFMs display interesting properties in their own right, including evidence of quantum tunnellingi^ 
Clearly there is a wide range of phenomena in AFM spintronics that is yet to be explored. 

VI. SPINTRONICS IN MOLECULAR SYSTEMS/PSEUDOSPINTRONICS 

In the previous two sections, new and interesting physics was revealed by considering CIT effects in materials other 
than FM (ferrimagnets and antiferromagnets) . That provokes the question of what other types of materials or systems 
might provide interesting manifestations of CIT. 

One type of system in which CIT effects may occur in interesting ways is in molecular systems, or those with reduced 
dimensionality. So far experiments have focused on magnetoresistance effects for systems with ferromagnetic leads 
sandwiching monolayers of molecules^ or carbon nanotubes;^ for example. The non-equilibrium Green's function is 
well suited to calculating spin-dependent transport for such systems<^ A next step would be to consider the action 
of current on magnetic molecules to see if current-induced switching of the molecular spin is possible, and to consider 
experimental signatures of such a switching event i^i^ Experimental challenges are abundant, and certainly outside 
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of the expertise of these authors, however first-principles calculations may be useful in aiding the identification of the 
optimum experimental choice for molecule and lead material. There can however be difficulties in applying density 
functional theory (strictly speaking a ground state mean-field theory) directly to molecular transport (a non-ground 
state system in which electron correlations may be important) It is nevertheless likely that the current-induced 

torques in such systems will be qualitatively different than in conventional FMs. In FMs, the properties of current- 
induced torques rely critically on the dimensionality (the spin-current becomes aligned to the local magnetization 
only after averaging over all incoming electron velocity directions) , and these molecular systems are effectively 1-d or 
0-d. 

Finally we mention that current-induced torques are closely related to physics that occurs in other kinds of systems. 
For example, the supercurrent that flows through a superconductor in a circuit with a normal metal source and drain 
can be thought of as being driven by current-induced changes in the equation of motion for the order parameter. In 
this case, total particle number (a scalar) rather than total spin (a vector) is conserved. The case of circuits containing 
superconductors is therefore closely analogous to the case of XY easy-plane ferromagnets for which only the z-spin is 
conserved. Bilayer quantum Hall systems near filling factor v ~ 1, have an exciton condensate^ ground state with 
spontaneous phase coherence in the two layers. The conserved quantity in this case is the difference between the 
numbers of particles in the two layers. The anomalous transport properties of these exciton superfluids are closely 
related^ to the anomalous transport properties of magnetic systems that we have discussed in this review. We suspect 
that many other examples of current-driven order parameters will be discovered and exploited in the future. So far, 
metallic magnetism has provided the most phenomenologically rich and most extensively explored example of this 
type of physics. The lessons learned from this still developing body of research may have implications beyond the 
realm of magnetism. 

VII. CONCLUSION 

In this review, we have described a microscopic theory of current-induced torques which identifles the interaction 
between misaligned spins of non-equilibrium quasiparticles and the magnetic condensate as the source of the torque. 
This perspective suggests that the phenomena of current-driven order parameters is more general than the spin- 
transfer idea. It suggests that in magnetic systems the phenomena is not limited to ferromagnets, or to systems in 
which total spin is conserved. We have applied this picture to consider spintronics in antiferromagnets, and have found 
a number of qualitative differences in the physics of CITs compared to the ferromagnetic case. There are a number 
of new areas that have not yet been fully explored with this approach - among them systems with strong spin-orbit 
coupling, different types of antiferromagnetic systems, and molecular systems. It is our hope that in exploring these 
novel systems, the key, fundamental aspects of current-induced torques can be illuminated in their most general form. 
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The functional derivative SE[rh]/{Srh) is defined by the energy variation SE = J dx[(5£'/((5m(x))] ■ 5m(x) of the energy 
functional. Writing the energy in terms of an energy density S{m), i.e., E[rh] = J dyi£{rh{x)), we have that 

SE[m] _ d£ d d£ 



5rh dm dx~f d{drh/dx^) 



where a sum over the repeated index 7 £ {x, y, z} is again implied. See also the introductory article for a discussion of 
functional derivatives. 



